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A simple and economical procedure for accurate determinations of toughness and lifetime parameters of 
ceramics is described. Indentation flaws are introduced into strength test pieces, which are then taken to failure 
under specified stressing and environmental conditions. By controlling the size of the critical flaw, via the 
contact load, material characteristics can be represented universally on "master maps" without the need for 
statistical considerations. 

This paper surveys both the theoretical background and the experimental methodology associated with the 
scheme. The theory is developed for "point" flaws for dynamic and static fatigue, incorporating load explicitly 
into the analysis. A vital element of the fracture mechanics is the role played by residual contact stresses in 
driving the cracks to failure. Experimental data on a range of Vickers-indented glasses and ceramics are included 
to illustrate the power of the method as a means of graphic materials evaluation. It is demonstrated that basic 
fracture mechanics parameters can be measured directly from the slopes, intercepts and plateaus on the master 
maps, and that these parameters are consistent, within experimental error, with macroscopic crack growth laws. 
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parameters. Chief among these parameters are the frac- 

Introouction ture toughness, K c , and the crack velocity exponent, n , 

The increasing use of glasses and ceramics as struc- which respectively characterize the equilibrium and ki- 

tural materials has prompted the development of new netic crack growth responses. In the context of brittle 

and accurate techniques for evaluating intrinsic fracture design it is essential to. achieve an adequate level of 

precision in such parameter evaluations. This is particu- 
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problem in any truly independent way. Evaluation of ma- 
terial parameters becomes a mere exercise in statistical 
data manipulation, with little or no physical insight into 
the nature of the critical flaws responsible for failure 
[1-2]. 1 This probabilistic approach makes it difficult to 
assess the relative merits of different materials from the 
standpoint of intrinsic properties alone. 

A controlled-flaw technique which effectively elimi- 
nates the statistical component from strength testing has 
been described in a series of recent articles [3-12]. A 
single dominant flaw of predetermined size and geome- 
try is introduced into the prospective tensile surface of 
each specimen using a standard diamond indenter. The 
specimens are then stressed to failure in the usual way. 
With the indentation and flexure testing conditions held 
fixed, any variations in the strength behavior can be 
taken as direct reflections of the intrinsic material re- 
sponse. The only need for statistical treatments then 
resides in the trivial accountability of random scatter in 
the data. Quite apart from the ensuing improvements in 
data reproducibility, the indentation procedure confers 
several advantages in strength analysis: (i) greater speci- 
men economy, (ii) because the location of the critical 
flaw is predetermined, closer observation of the fracture 
mechanics to failure, and (iii) a reasonable simulation of 
the damage processes that are responsible for a great 
many brittle failures [13-15]. One apparent complication 
attending the technique is the existence of a strong re- 
sidual contact field about the elastic/plastic deformation 
zone, necessitating the incorporation of additional terms 
in the governing stress intensity factor. However, 
closed-form solutions of the fracture mechanics formu- 
lations are now available for both equilibrium [4] and 
kinetic [16] conditions of failure; analytical deter- 
minations of toughness and fatigue parameters from the 
strength data may accordingly be made in as straight- 
forward a manner as for "Griffith" flaws without the 
residual stress term. 

The capacity to control the scale of the critical flaw 
via the indentation load is a potent tool in the in- 
vestigation of material fracture properties. The load ac- 
tually replaces initial crack size as a variable in the frac- 
ture equations, thereby eliminating the need for onerous 
measurements of crack dimensions (although some ob- 
servations of crack growth are useful for confirming the 
validity of the theory) [15]. Size effects in the micro- 
mechanics may then be studied systematically: im- 
portant changes in the nature of low-load contact flaws 
have been thus revealed on reducing the crack size to 



the scale of the deformation zone [17] or of the micro- 
structure [18]. Systematic variations in the load de- 
pendence of indentation-strength characteristics can 
also be used to evaluate pre-existing stress states in brit- 
tle materials, e.g., in tempered glass [19]. Again, some 
materials may produce ill-defined indentation patterns 
outside certain ranges of flaw size, or be restricted in 
specimen dimensions, in which case the geometrical re- 
quirements of standard strength-testing procedures may 
make it impossible to operate at a single contact load. 
The theoretical analysis allows one to compensate for 
any such changes in the working contact conditions, 
effectively reducing all data to an "equivalent" load. 
This paper illustrates a procedure for representing the 
intrinsic strength properties of brittle materials on an 
indentation "master map." A suitable "normalization" 
scheme incorporating indentation load into the plotting 
coordinates allows for the reduction of all inert and 
fatigue strength data on to "universal" curves for the 
various test materials. In this sense the scheme is remi- 
niscent of that developed earlier by Mould and South- 
wick [20], except that their use of relatively ill-defined 
abrasion flaws necessitated a totally empirical approach 
in the data reduction. On our master map the position of 
a given curve may be taken as a graphic indicator of the 
intrinsic toughness and fatigue susceptibility. Quan- 
titative determinations may accordingly be made of K c 
and n without recourse to statistically based theories of 
strength. 



Background Theory 

Stress Intensity Factor for Indentation Cracks 

The starting point in the analysis is the stress intensity 
factor for an indentation crack of characteristic dimen- 
sion c produced at peak contact load P and subjected to 
subsequent applied tensile stress 0%. For "point" flaws 
produced by axially loaded indenters the general form 
of this stress intensity factor is [4] 



K = X P/c 1/2 +x{iar 3 c 



1/2 



(1) 



where x an d $ are dimensionless parameters. The sec- 
ond term in eq (1) is the familiar contribution from the 
ap*plied field; ii depends only on crack geometry, here 
assumed to be essentially "penny-like" [21]. The first 
term is the contribution from the residual contact field; 
for materials which deform irreversibly by a constant 
volume process 



'Figures in brackets indicate literature references at the end of this 
paper. 



X=%E/H) 



1/2 



(2) 
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approximately [22], where E is Young's modulus, H is 
hardness and £ is a numerical constant. 

In the event of any pre-existent stress acting on the 
crack a third term would have to be included in eq (1) 
[4,9]. Other than to note that this potential complication 
needs to be heeded when preparing the surfaces of test 
specimens we shall consider it no further in our mathe- 
matical derivations. 



Equilibrium Solutions: Inert Strengths 

Equilibrium conditions of crack growth are closely 
realized experimentally by testing in an inert environ- 
ment. In terms of fracture mechanics notation the crite- 
rion for equilibrium is that K=K Q ; if dK/dc<0, the 
equilibrium is stable, \idK/dc >0, it is unstable. Now it 
is evident from eq (1) that K for given values of P and cr a 
passes through a minimum in its functional dependence 
on c; thus at subcritical configurations K(min)<.K c 
there is a stable and an unstable equilibrium, to the left 
and to the right of the minimum, respectively [16]. In an 
inert strength test, cr a is increased steadily until these 
two equilibria merge at dK/dc=0, which defines the 
critical variables 



cr m = 3K c /4xijc m l/2 



(3a) 
(3b) 



at which crack growth proceeds without limit. We may 
note that any relaxation of the residual stress field, as 
reflected in a reduction in x (or, more specifically, in £ in 
eq (2)), will cause <r m to expand and c m thence to con- 
tract. 

It can be shown that the ideal indentation crack is in 
a state of equilibrium immediately after completion of 
the contact cycle [22]. The size of this crack is found by 
setting o- a = 0, K=K C in eq (1); 



c,={ X P/K Q f 



(4) 



From eq (3b) we have c =0.40c m . On subsequently ap- 
plying the tensile stress, the crack extends stably from c 
to c m , whence spontaneous failure ensues at cr^o^n [4]. 
In reality, deviations from this ideal behavior are ob- 
served; relaxation effects can cause c m to contract, as 
already mentioned, and subcritical, moisture-assisted 
crack extension within the residual contact field can 
cause Co to expand, to Co say. Nevertheless, unless the 
condition Co<c m is violated, some precursor crack 
growth will still precede failure, in which case o- m re- 
mains a measure of the inert strength. 

Equation (3) may then be conveniently rearranged to 



eliminate all terms in crack size, and then combined with 
eq (2) to yield 



a- m P w ={yA^){\/^) Ui [{H/E) m K,] 



s\/lV 14/3 



(5) 



This expression conveniently relates the test variables 
on the left side to the material properties, primarily the 
toughness, on the right side. We emphasize once more 
that this formulation is contingent on the absence of all 
spurious pre-present stresses. 



Kinetic Solutions: Dynamic Fatigue 

When cracks are exposed to moisture or other inter- 
active environmental species, extension can occur in the 
subcritical region, K<K C . The major characteristic of 
this kind of extension is its rate dependence, which is in 
turn highly sensitive to the crack driving force. The 
basic equation of kinetic fracture accordingly takes the 
form of a crack velocity v (K). In the interest of obtain- 
ing closed-form solutions to the ensuing fracture me- 
chanics relations we choose the empirical power -law 
function [23] 



v=v (K/K c y 



(6) 



where v and n are material/environment parameters. 
Materials with lower values of n are said to be more 
"susceptible" to kinetic crack growth effects. 

The most practical loading arrangement for the sys- 
tematic study of rate effects in strength properties is that 
of "dynamic fatigue," in which the time differential of 
stress is held fixed up to the point of failure, i.e., 
&a=o-Jt =const. We may thus combine eqs (1) and (6) 
to obtain a differential equation for this stressing config- 
uration, 



dc/dt =v [xP/K c c V2 + xl>&*c U2 t /K c ]" 



(7) 



This equation has to be solved at given P and 6- a for the 
time to take the crack from its initial configuration, 
K =K(ci), to its final configuration, K=K C , at which 
point the stress level defines the dynamic fatigue 
strength, o- a = o- f [16]; 



where 



<r f =(^a) w+!) 

«' = 3«/4+l/2 
X' = (27r«') 1/2 o- in "' Cm /vo. 



(8) 

(9a) 
(9b) 
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The solution in eq (8) is identical in form to that for 
"Griffith" flaws (x=0) [23]. However, the slopes and 
intercepts from a linear plot of log cr f against log cr a are 
very different in the two instances. In the present case 
(v^O), n ' and X' may be regarded as "apparent" fatigue 
parameters, in the sense that transformation equations 
are required to convert these to "true" crack velocity 
exponent and coefficient terms. Thus, eq (9a) may be 
inverted to obtain n directly from n ', and eq (9b) simi- 
larly (in conjunction with measured values of <r m and c m ) 
to obtain v from X'. It is again seen that initial crack size 
does not enter the results, as long as the condition Co<c m 
remains operative [9]. 

Implicit in the derivation of eq (8) is the usual assump- 
tion that the prospective test surfaces are free of spuri- 
ous stresses. The introduction of such stresses leads to 
nonlinearities in the dynamic fatigue plotting scheme, 
thereby destroying the basis for the above analysis 
[9,10]. 

It is convenient at this point to incorporate the inden- 
tation load as a working test variable into the dynamic 
fatigue relations. Whereas n ' in eq (9a) is independent of 
all test variables, X' in eq (9b) can be expressed as an 
explicit function of P via the quantities cr m and c m in eq 
(3). In this way we may write 



Kinetic Solutions: Static Fatigue 

Of more practical interest from a design standpoint is 
the issue of component lifetime under fixed stress rather 
than stress rate . Ideally, it would seem desirable to for- 
mulate a universal static fatigue relation in direct anal- 
ogy to eq (12) retaining, as far as possible, the same 
adjustable parameters. Lifetime predictions could then 
be made from dynamic fatigue data alone, without hav- 
ing to resort to delayed failure experiments. This formu- 
lation may be achieved in two steps. First, eliminate 
stressing rate in favor of time to failure, 6- a =cr f /? f . This 
step introduces the lifetime concept without yet altering 
the status of eq (12) as a dynamic fatigue relation. Then, 
convert to equivalent static fatigue variables by re- 
placing o- f with o- A , i-e., the level of the invariant applied 
stress, and t t with (n' + l)t f [16]. The resulting static fa- 
tigue relation is 



U/P 2n = X' P /{n , + \){cr A P ui r 



(13) 



X'=X' P /P ( "'- 2)/3 



(10) 



We reiterate here, at the risk of laboring the point, 
that the variables P, cr A and t f in eq (13) relate to pro- 
spective static fatigue conditions, whereas the parame- 
ters n ' and X P are adjustables, as defined by eqs (9) and 
(10), to be determined from dynamic fatigue data. 



where X' P is a modified intercept term, totally indepen- 
dent of P, given by 

X' P = (2tt« ') 1/2 (3i(: c /4i|/)"'(ii:c/4x) ( "' -2)/3 /vo. (1 1) 

Equation (10) tells us that fatigue data obtained on one 
material but using different indentation loads will fall on 
different straight lines, mutually translated but without 
change of slope. Now by inserting eq (10) into eq (8) we 
may appropriately modify the dynamic fatigue relation, 
thus 



cr i P m = {X , P &J>) 



l/(n'+l) 



(12) 



such that by plotting log (o- f P l/3 ) against log (&*P) all 
data should fall on to a universal fatigue curve. This plot 
would, of course, cut off at a limiting level on the ordi- 
nate corresponding to the inert strength plateau defined 
in eq (5). The procedure for evaluating crack velocity 
parameters from the slopes and intercepts of such repre- 
sentations is the same as before, but with eq (10) serving 
as an intermediary to eq (9). 



Experimental 
Materials Selection and Preparation 

The materials in this study were chosen in accordance 
with two major criteria: first, they should cover a range 
of toughness and crack velocity characteristics, as deter- 
mined by independent fracture techniques; second, they 
should be of some technical importance. Table 1 lists 
these materials and their pertinent properties. 

All specimens were prepared in the usual manner for 
strength testing. However, particular attention was paid 
to surface preparation, bearing in mind our repeated 
assertion that pre-existing stress states can greatly influ- 
ence the interpretation of strength data. The glass speci- 
mens were therefore annealed [19] and the ceramics 
surface-polished to a mirror finish with diamond paste 
[10] to ensure removal of any such stresses. 



Indentation and Strength Testing Procedure 

All specimens were routinely indented centrally 
along their length using a Vickers diamond pyramid 
indenter to produce dominant flaws for the subsequent 
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Table 1. Materials used in this study. 





Independent parameters 


Indentation 


parameters 






E 


H 


K c 


n 


K< 


n 


log v 


Materials 


















(sources footnoted) 


GPa 


GPa 


MPam ! 






MPam ! 




ms" 1 


Soda-Lime Glass" 


70 


6.6 


0.74* 


16- 


19* 


0.97 


18 


-1.6 


Borosilicate Glass b 


89 


6.5 


0.77» 


31- 


37* 


1.2 


36 


1.6 


Fused Silica 


72 


7.6 


0.81* 


38* 




1.2 


44 


2.2 


Synroc d 


190 


10.3 


1.9 






1.8 


35 


0.2 


P.Z.T." 


88 


3.1 


0.87 






1.0 


43 


-0.5 


Alumina f 


400 


16 


4.4 


46* 




3.8 


59 


1.7 


Silicon Carbide 5 


435 


24 


4.1* 


118" 




3.7 


222 


8.4 


Glass Ceramic 11 


108 


8.4 


2.5* 


63,* 


84* 


2.2 


117 


5.0 



'Determinations by other workers. (See references, below). 

a Schott-Ruhrglas GMBH [11,31,32]. 

b Schott-Ruhrglas GMBH [11,31,32]. 

c Schott-Ruhrglas GMBH [31,33], 

d Synroc B, Australian Atomic Energy Research Establishment [34], 



c Lead Zircon Titanate, Plessey Australia. 

r F99, Friedrichsfeld GMBH [35]. 

BNC203, Norton Co. [7,36], 

h Pyroceram C9606, Corning Glass Co. [7,10,37,38]. 



failure tests. The Vickers geometry was chosen both for 
its proven capacity to produce well-defined radial crack 
patterns and for its general availability in hardness test- 
ing facilities. The glasses were indented at several loads, 
ranging from 0.05 to 100 N, whereas the ceramics were 
each indented at single loads, 10, 20, or 100 N. In all 
cases the radial cracks extended well beyond the central 
hardness impression, but never to a length in excess of 
one tenth the specimen thickness. 

The indented specimens were then broken in four- 
point flexure [24] in a universal testing machine at con- 
stant crosshead speed. Care was taken to center the 
indentation on the tension side, with one set of radial 
cracks aligned normal to the long axis. The breaking 
loads were recorded using conventional strain gage and 
piezoelectric load cells [10], and the corresponding rup- 
ture stresses thence evaluated from simple beam theory. 
Inert strengths, cr m , were measured in dry nitrogen or 
argon or silicone oil environments, with the crosshead 
running at its maximum speed. Dynamic fatigue 
strengths, cr f , were measured in distilled water over the 
allowable range of crosshead speeds. At least six speci- 
mens were broken in each strength evaluation, from 
which means and standard deviations were computed. 



Measurement of Critical Crack Dimensions 

For the purpose of confirming the necessary condi- 
tion that the initial crack size Co should never exceed the 
instability value c m for equilibrium failure, and for veri- 
fying certain aspects of the fatigue solutions presented 
earlier, an optical examination of representative critical 



indentations is recommended. The technique used here 
was to place three indentations instead of one on a given 
test surface, and then take the specimen to failure under 
inert conditions [10]. On the understanding that all three 
indentations must have had nearly identical growth his- 
tories, the procedure leaves two "dummies" in the bro- 
ken test piece from which to measure the required crack 
dimensions. The Vickers geometry proves particularly 
useful in this technique, for while the set of radial cracks 
perpendicular to the tensile direction provides a mea- 
sure of c m , the set parallel to this same direction remains 
free of external stress and hence provides a measure of 
Co. 

In all materials studied in this work some precursor 
crack growth was indeed found to occur prior to failure. 



Results 
Inert Strengths and Toughness 

In this section we begin by examining the dependence 
of inert strength on indentation load for the three glasses 
studied. With this dependence established, we then in- 
vestigate how the inert strength data may be reduced to 
a composite toughness parameter for all of the test mate- 
rials. 

Figure 1 accordingly shows <x m as a function of/ 3 for 
the glasses. The straight lines are best fits of slope — 1/3 
in logarithmic coordinates, as per eq (5). This same de- 
pendence has been confirmed elsewhere for several 
other brittle materials [7,18,25,26]. 

Values of the composite parameter cr m P m are thus 
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Figure 1 — Inert strength as function of indentation load for the silicate 
glasses. (Data courtesy of T. P. Dabbs). 

evaluated for each of the glasses and ceramics, and are 
plotted as a function of (H/E) m K c (from table 1) in 
figure 2. The straight line is a fit of logarithmic slope 4/3 
in accordance with eq (5), using a "calibration" value 
(3/4il/)(l/4§) 1/3 = 2.02 from an earlier, more comprehen- 
sive study [7]. The trends in figure 2 appear to be in 
reasonable accord with prediction, although some devi- 
ations are evident, particularly for the fused silica and 
borosilicate glasses. Estimates of the "indentation 
1000 
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Figure 2— Inert strength parameter cr m P m as function of toughness 
parameter {H/E) m K Q for the glasses and ceramics. 



toughness" obtained directly from cr m P yi by inverting 
eq (5) are included in table 1 for comparison with the 
independently determined values. 



Dynamic Fatigue and Crack Velocity Parameters 

We consider now the dynamic fatigue responses, 
again beginning with the glasses to examine the func- 
tional influence of contact load, and outline the pro- 
cedure for determining the exponent and coefficient in 
the crack velocity function. 

Figure 3 shows these responses for the glass com- 
positions in water. The straight lines drawn through 
individual sets of data at fixed P are best fits to eq (8), 
regressed for each glass on all the data consistent with 
the intercept relation eq (10). Thus we obtain families of 
lines of constant slope, with systematic displacements to 
lower strength levels with increasing load. Analogous 
plots are shown in figure 4 for the five ceramics in the 
same water environment, but now for a single load in- 
each case. The inert strength limits are included in all 
plots as a reference baseline for assessing the degrees of 
fatigue. 

From the regressed slopes and intercepts we obtain 
values of the apparent fatigue parameters n ' and X' in eq 
(8). Inversion of eq (9) (together with the inert strength 
data) then allows us to evaluate the true crack velocity 
parameters, n and v . These evaluations are summarized 
in table 1; comparisons may be made in this tabulation 
with independent measurements of the crack velocity 
exponent. 



Master Maps 

We have set the base for determining universal frac- 
ture curves for the materials studied, and thence to con- 
struct master maps. We do this for dynamic and static 
fatigue conditions in turn. 

The presentation of the dynamic fatigue results on a 
single master map requires conversion of all data to 
appropriate load-adjusted variables cr { P ul and &J> in eq 
(12). Figure 5, an appropriate composite of all data thus 
converted from figures 3 and 4 (but with error bars 
omitted for clarity), is such a map. Each material is now 
conveniently represented by a universal curve, indepen- 
dent of the contact loads used to obtain the data. The 
curves plotted in this diagram represent numerical solu- 
tions of the basic fatigue differential equation, eq (7), 
obtained for the ranges of P and 6- a covered experi- 
mentally for each material, using the inert and kinetic 
parameters already determined along with the measured 
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initial crack sizes [10]. The fact that the curves regen- 
erated in this way are effectively coincident with the 
data is, of course, no real surprise, since the regression 
analyses used in the parameter evaluations were per- 
formed in accordance with the solutions of the differ- 
ential equation in the first place. An exercise of this kind 
nevertheless serves two useful purposes: (i) to confirm 
that the solutions referred to, which are of closed form, 
are indeed reasonably reliable, and (ii) to show how 
closely the curves remain linear in the fatigue region, 
and then plateau out at the inert strength levels, cr m P Ui 
(fig. 2). 

The equivalent construction for static fatigue is ob- 
tained from the constant stressing rate results using the 
rationale described earlier in the derivation of eq (13). 
Thus we generate the plots shown in figure 6 directly 
from the best-fit values of n ' and X' (or more strictly, via 




10- 2 



10° 10 2 10 4 

Stress Rate, o a /M Pas -1 



10 6 



10 a 



Figure 3— Dynamic fatigue responses of glasses indented at different 
loads. The hatched bands indicate inert strength levels. (Data cour- 
tesy of T. P. Dabbs). 



eq (10), X'p) determined by the data regressions shown in 
figures 3 and 4. Cutoff levels on the abscissa again corre- 
spond to inert strength limits. Because the construction 
in figure 6 is not obtained this time from regenerated 
solutions of the basic differential equation, we are unable 
to plot the curved transition between the fatigue and 
inert regions; however, the abruptness of the corre- 
sponding crossover points in figure 5 suggests that we 
may reasonably ignore any such curvature in the life- 
time maps. 

Discussion 

Quantitative Evaluation of Fracture Parameters 

The scheme presented here for reducing fatigue data 
to universal curves for any specified material/ 
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Figure 4-— Dynamic fatigue responses of ceramics indented at single 
loads. Hatched bands indicate inert strength levels. 
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Figure 5— Dynamic fatigue master 
map. Error bars omitted from 
data points for clarity. See table 
1 for key to materials. 
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environment system, and thence for constructing master 
maps to facilitate comparisons between these curves, 
provides an attractive route to simple, accurate, and 
economical evaluation of fracture parameters for de- 
sign. In the following subsections we discuss how these 
constructions may be used as a quantitative tool for 
parameter determinations in different regions of the 
curves. 

Inert Strength Levels. The position of the inert 
strength cutoff level, cr m P m , may be taken as an indi- 
cator of material toughness, K c . Intrinsically tougher 
materials will therefore exhibit cutoffs further toward 
the top of a dynamic fatigue map (fig. 5) and toward the 
right of a static fatigue map (Fig. 6). 

It should be emphasized that the correspondence im- 
plied here is not exact. To clarify this point we may 
invert eq (5) to obtain an explicit expression for tough- 
ness, 



K c = (256^/27) m (E/H) m (cr m P ,/3 ) 



l/3\3/4 



(14) 



Thus, K c depends on the elastic/plastic term E/H as 
well as on cr ra P 1/3 . On the other hand, since E/H varies 
only between 10 and 25 over the range of materials listed 
in table 1, the use of an invariant, representative mean 
value <{E/H) Ui > = 1.50 in eq (5) would lead to errors 
of no more than 10%. Another potential source of dis- 



crepancy lies in the implicit assumption that geometrical 
similarity is preserved in the indentation pattern from 
material to material, as reflected in the constancy of the 
parameters £ and ty. We have already pointed out that 
relaxation effects in the residual contact field can lead to 
reductions in the £ term. Systematically low values of £ 
will also be manifest in materials which deform by other 
than a constant- volume process or exhibit plastic pile up 
at the impression edges [22]. Fused silica and boro- 
silicate glass, which tend to deform by densification [27], 
fall into this category, thereby explaining the tendency 
for the data points representing these two materials to lie 
above the general trend in figure 2. Finally, it has been 
taken as given that the radial crack patterns are always 
well defined, and in the .materials used here this gener- 
ally has been found to be so. But in materials where the 
microstructure is comparable in scale with the inden- 
tation event, the symmetry of the crack pattern can 
become severely disrupted [8,28], with consequent vari- 
ations in both § and i|/. 

It may be argued that the "effective" toughness reck- 
oned from the cutoff position on a master map, while 
perhaps not an accurate measure of its macroscopically 
determined counterpart, nevertheless may more closely 
characterize the response of "natural" flaws. This is 
certainly likely to be so where the strength-controlling 
flaw in a component is created by a surface contact 
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event, as in sharp-particle impact or in a machining 
operation. In this sense the master map approach might 
well be expected to serve more appropriately as a source 
of design parameters than the more conventional meth- 
ods involving large-scale fracture specimens. 

Fatigue Curve Slopes. We have noted from eqs (12) 
and (13) that the slope of a universal fatigue curve is a 
measure of the intrinsic susceptibility to slow crack 
growth. Thus materials with lower values of the crack 
velocity exponent n , and hence of n ', eq (9a), will have 
greater slopes on dynamic fatigue master maps and, con- 
versely, lower (negative) slopes on static fatigue maps. 

As with the toughness, certain caution needs to be 
exercised when using master map data to determine n 
values. This is because in applying the inverted form of 
eq (9a), 



«= 4b 73 -.2/3, (15) 

it is implicit that certain necessary conditions are met. 
The most important of these is the proviso Co<c m , which 
we have considered at some length in this work. It is 
interesting to note that if this proviso is satisfied even the 
"anomalous" glasses which deform by non-volume- 
conserving processes may be analyzed in terms of eq 
(15); the fatigue properties are not sensitive to the origin 
of the residual contact field, as long as this field is of 
sufficient intensity to generate some precursor crack 
growth [11]. If such a precursor stage were not to be 
evident in the failure mechanics the "apparent" term n ' 
would tend closer in value to the true n [5,9,12,13]. A 
second condition that needs to be met is that the flaws 
should indeed be produced in axial loading; other inden- 
tation loading systems, e.g., linear translation, give rise 
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to flaws which are governed by a transformation equa- 
tion with coefficients significantly different from those 
in eq (15) [12,16]. 

It is seen in table 1 that the exponents obtained from 
this study agree well with the independent deter- 
minations for the glasses, but not for the ceramics. The 
relatively good agreement in the case of the glasses is 
attributable in part to the "model" behavior of this class 
of materials: transparency, isotropy, absence of micro- 
structural complication, and ease in specimen prepara- 
tion are factors which contribute to this behavior. Also, 
the n values of the glasses are comparatively low, so 
fatigue effects show up more strongly. This last point, 
coupled with a growing realization that conventional 
testing techniques used to obtain macroscopic velocity 
data are themselves subject to uncertainty (particularly 
the double torsion specimen [29]), could account for the 
discrepancies evident in the data for the ceramics. 

Fatigue Curve Intercepts. The intercept terms in the 
master map representations do not have such a simple 
interpretation in terms of basic fracture parameters. 
This is clear from eq (11); X P is a function of several 
quantities. Given the fatigue slope and inert strength 
evaluations as described in the two previous subsec- 
tions, along with a direct measurement of the critical flaw 
size c m , \ P effectively determines the crack velocity co- 
efficient v . Due to the compounding of errors (particu- 
larly from the n ' exponent), determinations of this kind 
are subject to gross uncertainty. There accordingly 
seems to be little value in trying to retain v as a design 
parameter, particularly since the X' terms, which can 
usually be determined to within 15% from dynamic fa- 
tigue data, may be used directly in lifetime formulae. In 
studies of the basic physics and chemistry of crack 
growth, of course, v remains a useful coefficient for 
scientific analysis. 



Practical Implications of Master Maps 

The major appeal of the master map construction ad- 
vocated here lies in the provision of a graphic indicator 
of the intrinsic toughness and fatigue properties of brit- 
tle materials. Each material is represented by a universal 
curve, the relative position of which determines the 
merit of that material for structural applications. The 
marked superiority of such materials as silicon carbide 
and alumina become vividly apparent in the maps of 
figures 5 and 6. Useful distinctions may also be made 
between materials which cross over within the data 
range, e.g., soda-lime glass and PZT. On the basis of 
straight inert-strength testing we might reckon the first 
of these as the stronger material, whereas for applica- 



tions involving sustained stresses it is the second which 
would tend to the larger lifetimes. Such crossovers 
would not be so obvious from the raw fracture mechan- 
ics parameters. It will be appreciated that this kind of 
intercomparison is made on the basis of "equivalent" 
flaw sizes: in this respect the indentation method, 
through its control over the flaw severity via the con- 
tact load, is unique in its capacity for reducing strength 
data to a common denominator. 

In arguing the merits of this approach we do not mean 
to imply that it is only the intrinsic fracture properties 
which play an important role in the determination of 
component strengths and lifetimes; the effective sizes of 
the naturally occurring flaws which ultimately cause 
failure must also be known. Our procedure, by intro- 
ducing flaws greater in severity than any of these natural 
flaws, automatically excludes information concerning 
the latter from the data. What our scheme effectively 
allows us to do is to determine the intrinsic parameters 
in a truly independent manner. All necessary extrinsic 
flaw parameters should be obtainable from straight- 
forward inert strength tests (run at a single stressing 
rate), in the form of the usual statistical distribution 
functions. Lifetime predictions for as-prepared com- 
ponents could then be made without ever having to 
accumulate vast quantities of fatigue data [2]. In adop- 
ting this strategy one needs to keep in mind the strong 
influence that any persisting residual stress concen- 
trations associated with the original initiation processes 
(in our case the elastic/plastic deformation) might exert 
on the subsequent flaw evolution. In the absence of in- 
formation as to this aspect of flaw characterization steps 
should be taken to design conservatively, on the basis of 
"worst-case" configurations wherever possible. This 
last point is dealt with in greater detail in reference [25]. 

It has been indicated at several points that the exis- 
tence of any spurious stresses incurred during the me- 
chanical, chemical, or thermal history of a material 
would necessitate a third contribution to the starting 
stress intensity factor in eq (1), with consequent devi- 
ations from the currently determined toughness and fa- 
tigue relations. The fact that such deviations were not 
observed in the materials studied here may be taken as 
evidence that this potential complication has been suc- 
cessfully avoided. Again, it may be well to emphasize 
that it may not be so simple to confirm the elimination of 
spurious stresses from surfaces whose strengths are con- 
trolled by natural flaws, particularly in materials with 
typically wide flaw distributions; nor, of course, may we 
wish to eliminate them, bearing in mind that these 
stresses are most often compressive. 

Finally, a comment may be made concerning the con- 
venience of indentation load as a variable for in- 
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vestigating fundamental flaw size effects. By system- 
atically reducing the load we can produce correspond- 
ing smaller flaws, thereby providing a link between 
macroscopic and microscopic crack behavior. Any 
change in the nature of the indentation flaw will then 
become evident as deviations from universal plots, 
much as just described in relation to the spurious stress 
influence. In this way it has been possible to demonstrate 
that indentation flaws in glasses undergo an abrupt tran- 
sition in properties below a threshold load (correspond- 
ing to a flaw size ;= jam): above this threshold the macro- 
scopically determined laws of crack growth remain 
perfectly valid, regardless of scale, provided that the 
residual contact term is duly accounted for [11]; below 
the threshold the universal curves no longer apply, and 
failure becomes dominated by initiation micro- 
mechanics [17,30]. The indentation technique should 
prove similarly useful for studying size effects in ce- 
ramics, particularly for polycrystalline materials with 
relatively coarse microstructures. 



Conclusions 

1) The indentation-flaw technique provides an attrac- 
tive route to the evaluation of intrinsic fracture parame- 
ters. Coupled with independent determinations of natu- 
ral flaw distributions, the approach offers the prospect 
of accurate lifetime predictions with optimum specimen 
economy. 

2) The control over the nature, shape, and above all, 
the size (via the contact load) of the indentation flaw 
allows for the derivation of a universal fracture formu- 
lation. Each material is represented by a single curve 
which incorporates the toughness and fatigue proper- 
ties. Composite plots of these curves produce master 
maps, affording a simple graphic format for materials 
comparisons. 

3) The inert strength cutoff on such a master map is 
a measure of effective material toughness. For "well- 
behaved" materials this effective toughness is consistent 
with macroscopically measured K c values. In cases 
where inconsistency is observed the toughness reck- 
oned from indentation data may provide a more reliable 
indication of the response of the typical natural flaw. 

4) The slope of the fatigue curve on a master map is 
a measure of the susceptibility of a material to subcritical 
crack growth. The crack velocity exponent determined 
from this slope is an apparent value, « ', which is con- 
verted to the true value, n , via a simple transformation 
equation. 



5) Deviations from universality on a master map indi- 
cate an extraneous influence in the fracture mechanics, 
e.g., spurious stress states, microstructure/crack inter- 
actions, and threshold size effects. 



The authors thank T. P. Dabbs for providing raw 
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